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1. Introduction. 

In this paper we will investigate the singular points of the following unstable free 
boundary problem: 

(1.1) Au = - X {u>o} inBi(O) 

where X{u>o} is the characteristic function of the set {u > 0}. 

This problem was first investigated by G.S Weiss and R. Monneau [13] ■ In [T3] . 
C 1 -regularity locally energy minimising and maximal solutions of (jl.ip is shown. 
There is also some discussion regarding the possibility of the existence of singular 
points, that is points x° G -E>i(0) such that u ^ C 1 ' 1 (B r (x )) for any r > 0. Such 
points are proved to be totally unstable [14] . 

Let us formally define singular points before we proceed. 

Definition 1.1. Let u be a solution to Then we define S[u), the set of 

singular points of u, according to 

S(u) = {x G Bi(0); u i C 1,1 (B r (x)) for any r > 0} . 

Furthermore we will denote by S n -2{u) the singular points of co-dimension 2: 

u(rjX + y) 



i(u) = < y G S(u); lim - 

r 3 -^o \\u(rjTi + y)\\ L 2 (Bl ( 0)) 



r 2 



j-jj for some Q G Q and rj 



11^-1-^11^(51(0)) 
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where Q is the matrix group of rotations of K™ . 
It was shown in [13] or [3J that if y G S(u) then 

r u(r,-x + y) 

hm — f— - G V 2 

if the right hand side is defined, here V2 is the set of homogeneous second order 
harmonic polynomials of degree 2. Since the only homogeneous second order har- 
monic polynomial, up to translations, rotations and multiplicative constants, in R 2 
is xl — x\ it follows that S n -2 singles out the singular points with co-dimension 2 
singularities. 

In [3] two of the authors showed rigorously that singular points exists, that is 
there exist a solution u to (II. ip such that S(u) 7^ 0. This investigation was followed 
by the authors in [2] and [3J where we investigated and provided a total classification 
of singular points in R 2 and R 3 respectively. 

In this paper we intend to prove that in R n the singular points of smallest co- 
dimension are locally contained in a C 1 — manifold of dimension n — 2 and that the 
free boundary T u , defined 

r u = {xGBi(0); u(x)=0}, 

consists of two C 1 manifolds of dimension n — 1 intersecting orthogonally at such 
singular points. 

Our main theorem is 

Theorem 1.2. Let u be a solution to 11. 1]) and assume that 

u(rjx) x n-l ~ x n 



IK^ x )IU 3 (Bi) |kn-l -^nlU 2 (fli) 
for some sequence Tj — > (In particular, G S n -2{u) ). Then 

/, o \ ,■ u ( r j*) x\_ x - x 2 n 
(1.3) hm — + = 

and /or eac/i r\ > i/iere exists an r n > swc/i £/ia£ 

5ns n ,(0)n|x; |> 2 <»7(a£_i + a£)J 

consists of two C 1 manifolds intersecting at right angles at the origin. 
Furthermore there is a constant ro (u) > such that the set 

w(rx + y) 



S n -2 = S y; u(y) = |Vu(y)| = and lim 



r^o ||u(rx + y)|| L2(Bl 



Q 71 — 2 — ~ — T11 — ~ f or some Q ^ Q 



(0)) 



\\K,-i -4 IU 3 (si (o)) 

is contained in a C 1 manifold of dimension (n — 2) in i? ro (0). 

We would like to place this result in a long tradition of regularity result for 
parametric non-linear PDE. In particular we may view the free boundary T u — 
{x G Bi(0); u(x) = 0} as a parametric surface with singular points in S(u). 

Some of the most famous result in this area are the results by Bombieri, De 
Giorgi, Giusti and Simmons ([6], [17]) that states that no minimal cones exists for 
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minimal surfaces in n < 8. We should also mention the result by B. White [18] 
where uniqueness of tangent cones for 2-dimensional minimal surfaces is proved. 
From our point of view White's proof is interesting in that he uses a Fourier series 
expansion in constructing comparison surfaces. However, we work in n— dimensions 
which means that our Fourier expansions are considerably more subtle and involved 
than those that appear in [18) . 

Singularities in parametric problems have appeared in other areas of mathematics 
as well and our results have some similarities to the theory for harmonic mappings 
f [16j for a good overview). One could also mention a certain similarity with the 
theory of singularities that arise for a-uniform measures [13] . 

Equation ([l.lj) also arises in several applications for instance in solid combustion 
(see the references in [H]), the composite membrane problem ([8], [7], [5], [15], [9], 
[TO]), climatology ([H]) and fluid dynamics ([I]). 

Our proof will be based on a dynamic systems approach where we project a 
solution "^f^ onto the harmonic second order polynomials, call this projection 
II(it, r, 0) (see Definition 13. 2\i . By a careful analysis of the PDE we will be able to 
estimate H(u,r, 0) — II(u,r/2,0). Close to a singular point we have that "^"^ ~ 
H(u, r, 0) + ^n(ii,r,o) where 

AZ n (u,r,0) = — X{II(u,r,0)>0} m R™ 

Zn{u,r,o)(0) = |V^n(«,r,o)(0)| = 
hmixi^oo — L ^pH — = 

n(^n(«,r,0)j 1; 0) = 0- 

If we disregard lower order terms we may consider the map J-"(n(u, r, 0)) = 
n(M,r/2,0) defined by 

T(Tl(u, r, 0)) = IL(u, r, 0) + U(Z n(u>r>0) , 1/2, 0). 

The blow-up is unique if linifc^oo J- k (Jl(u, r, 0)) exists. 

Since the harmonic second order polynomials form a finite dimensional space. 
The map J 7 is a map between finite dimensional vector spaces. The main difficulty 
is that T is highly non-linear and we need quite subtle estimates to characterise 
the map. On the positive side we may write down Tl(u, r, 0) explicitly, modulo 
lower order terms, by means of Theorem 13.51 bv Karp and Margulis [12] , The 
definition of J 7 involves a Fourier series expansion of — Xn(«,r,o) 011 the unit sphere. 
Our main effort will be to estimate the Fourier coefficients in this expansion when 
n(w, r, 0)/ sup Bl | IT | w x 2 n _ x — x 2 n . For further details on the idea of the proof we 
refer the reader to [3]. 

2. List of notation: 

(1) 5 will denote a vector in M n_2 , we will always assume that \5\ << 1. We 
also define 5 = Yn=i 

(2) Vi will denote the second order homogeneous polynomials. 

(3) S(u) and S n -2(u) are the singular set and the singular set of co-dimension 
2 respectively, defined in Definition ll.il 

(4) The mapping F is defined in equation (I4.20p . 

(5) II(it, r, x°) is defined in Definition 13.21 

(6) The average of u in f2 will be denoted (u)q. 

(7) By dA we mean an area element of the surface under considration. 
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(8) We will use Landau's 0(r) notation to indicate a term that is bounded 
from by Cr for a universal constant C. That is f(x) = 0{r) if and only if 
\f{x)\ < Cr for a universal constant C. Similarly, f(r) > 0(r) means that 
f(r) > Cr for some universal constant C > etc. 

(9) ps(x) = YJiZi $i x i + (1 - &) x n-i ~ x ™> in particular p (x) = x^-i ~ x\. 

(10) Z Ps is defined in flUJ) . 

(11) Q is the matrix-group of rotations of W 1 . 

(12) The functions B 4 (5), 5(5), C»(5) and C(<5) are defined in (j4T2|) . (j4T3l) . 
Proposition 14. 31 and the remark after that Proposition respectively. 



3. Background Material and General Strategy. 

In this section we will state some of the results of [3] and outline our strategy 
(which is similar to the strategy of [3]). 

Our starting observation is the following proposition (Proposition 5.1 in |14| ) 

Proposition 3.1. Let u be a solution of in i?i(0) and let us consider a point 
x° G S(u). Then 

u(r,x + x°) „ 

lim j. — -, V 3 „.„ ' E V 2 

r 3 -H) \\u{rjX + x )\\ L 2 iBim 

for each sequence Vj —¥ such that the limit exists. 

The proof is a fairly standard application of a monotonicity formula. 

If u is a solution to (|1.1[) then Au G L°° which directly implies that D 2 u 6 
BMO(B 1 / 2 (0)) which in particular implies, via the Sobolev inequality, that for 
x° g S(u)nB VJ (0) 

(34) u(rx + x °) - l(x - x°)(J 2 U ) j? „ (x o ) (x - x°) 

is locally bounded in L 2 and pre-compact. It will be convenient for some calcula- 
tions later to subtract a harmonic polynomial in (|3.4[) instead of the polynomial 
i(x — x°)(£> 2 u) Br ( x o) (x — x°). We make the following definition. 

Definition 3.2. By II(u,r, x°) we will denote the projection operator onto V2 de- 
fined as follows: H(u, r, x°) = r r p, where r r 6 M + and p 6 V2 satisfies sup Bl \p\ = 1 
as well as 



inf 



D l [ M(rx + X°) \ _ D 2 h 2 _ f n 2 f "( rx + x °) 



'Bi(O) 

VFe wiH o/ten write II(u, r) w/ien x° is either the origin or given by the context. 
By definition T r — sup Bl |II(it,r)| and p r = H(u,r)/r r . 

It is a simple consequence of the BMO estimate (|3.4[) that if x° G S(u) H-Bi /2(w) 
then (Proposition 3.7 in [3]) 

u(rx + x°) 



(3.5) 



n(t»,r,x°) 



< C a (sup \u\,n). 

C7i,°(Si) Si 



If x° G S(u) then 



(3.6) sup \u\ > cr 2 ln(l/r) 

B r (x°) 
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for < r < ro(u,x°) and some small c > 0. To be more precise it is known that 
(c.f. Lemma 5.1 in [3]). 

Lemma 3.3. Let u be a solution to hl.l)) in B\ such that sup Bl \u\ < M and 
u(0) = |Vu(0)| = 0. Then there exist po > and r$ > such that if 

(3.7) sup|II(u,r)| > — 

Bi PO 

for an r < r then 

sup |n(u, r/2) | > sup |n(u,r)| + %A 

Bi Bj 

where t]q is a universal constant. 

The Lemma is proved for n = 3 in [3] but the proof is the same in arbitrary 
dimension. 

This estimate together with (|3.5j) implies that u(- + x ) = H(u, r, x°)+a lower 
order perturbation. Using the pre-compactness in C l a (c.f. Equation §3.5\i ) of 

(3.8) " (r,x 2 +x0) -n(v 3 ,x°) 

r j 

for some sequence rj — > we may extract a sub-sequence, which we still denote by 
rj , such that 

/ u(rjX + x°) rr , n .\ 
r % \1* -% r i. x ) = Z pW 

for some function Z p . It is not difficult to see that Z p is the unique solution to 

^ z p = -X{ P (x)>o} in M" 
^(0) = |VZ p (0)|=0 

(3 - 9) lim, , ^ - 

llm |x|— >00 | x |3 U 

n(z P) i) = o 

where 



p(x) = lim 



ri-m ||n(«,r J) xP)||£2( Bl )' 
In order to show regularity for the free boundary near a singular point we would 
have to control the limit 



lim 



n( M ,r J; x ) 



■r^o \\IL(u,rj,yfi)\\ L ^ Bl) ' 
If one can show that the limit is unique then it follows that the blow-up 
lim (u(rx + x°)/r 2 - U(u, r, x )) = Z D 

is unique. 

The following result, Corollary 7.3 in [3], gives a quantitative measure on how 
the function Zn(u,r,o) controls the difference between II(u,r, 0) and II(u,r/2,0). 

Proposition 3.4. Let u solve in B\ C K™ and assume that sup Bl \u\ < M, 

it(0) — |Vu(0)| = 0, and that for some p < po and r < r$, 

sup |II(u, r)| > — . 
Bi P 
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The 



sup |II(tt, r/2) - n(tt, r) - n(Z n( „, r ), 1/2)1 < C(M, n, a)(sup |n(«, r)|)" Q 

Bi Si 



/or eac/i a < 1/4. 

In order to estimate sup Bl ( ) |II(u, r, 0) — r/2, 0)| we thus need to be able to 
calculate H(Zu( u ,r,o): 1/2, 0). We will do this with the help of the following theorem 
from [12] 

Theorem 3.5. Let a G L°°(M. n ) be homogeneous of zeroth order, that is o"(x) = 
cr(rx) /or aZZ r > 0. Assume that a has the Fourier series expansion 

oo 

o-(x) = y^Q t crj, 

i=0 

on the unit sphere, where o~i is a homogeneous harmonic polynomial of order i. 

Moreover assume that AZ = a and that Z(0) — |VZ(0)| = lim^oo Z(x)/|x| 3 = 
0. Then 



Z(x) = g(x)ln|x| + |x| 2 0(x 



where 



a 2 

— r-^°"2 



<)=£ 



i x 

0"i 



2 (n + »)(i-2) 'Mx| y ' 
Our strategy in the rest of the paper will be to use Theorem 13.51 to calculate 
(3.10) n(Z n(u , r , ), 1/2, 0) = - k^<r a (x) 

where o~% is the second order term in the Fourier series expansion 

oo 

-X{n(«,r,o)>o} = ^OiO-^x) on &Bi(0). 

i=0 

Using the expression f|3 . 10[) in Proposition 13.41 will give us enough information to 
deduce that the blow-up of u is unique at all points x° 6 S n -2(u). 

4. Estimates of the Projections. 

In order to estimate U(Z p ,l/2) we need to calculate O20"2 from Theorem 13.51 
That involves calculating the second order Fourier coefficients for — X{p 5 >o} on the 
unit sphere. To that end we choose nxj — |x| 2 for i = 1, n and XiXj for i j as 
a basis for the second order harmonic polynomials. 

We may choose coordinates so that 

II(u r 0) 
SUPS! l n (u,r,0)| 

where 6 = (61,62, ^-2) and 6 — J^i^i ^i- We also define the polynomial ps, for 
a given vector 8 G K™~ 2 in equation (|4.11|) . We will assume, for definiteness that 
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5 > (this is implicit in the definition of ps in equation (|4.11|) ). If 5 < then all 
the following arguments follows through with minor and trivial changes. 

It follows from symmetry (i.e. — X{p s >o] i s even and the x^'s are odd on the 
unit sphere) that the Fourier coefficient of XiXj is zero. 

Since we are only interested in points x° G S n -2{u) where 

n(u,r j ,x") 

llm n FFF7 nTT = P°' 

r 3 — K) sup Bl \U(u,r 3 ,x u )\ 

for some sequence rj — > 0, we may assume that |<5 << 1. 
We also denote by Bi(5) the following integral 



(4-12) Bi(S) = - / X {p s >o}xidA 

and by B(S) the following integral 

(4.13) B(5) = - [ X{ PS >o}dA. 

Jasi(o) 

Here dA is the surface element. It follows that the Fourier coefficient of nx\ — |x 
of X{ Ps >o} is 

I, 2 i ,i, WS) - B(5j) . 

\\nxf - \\ L 2(d Bx (0)) 

Using that H(Z P5 , 1) = by definition and Theorem 13.51 we may deduce that 

n 

(4.14) n(Z„„l/2) = -K (n 2 Bi(S) - nB(S)) xj, 
where 

K ln(2) 



lL~(Bi(0))7 ' 



(n + 2)\\nx 2 ~ |x| 2 ||i2 ( a Bl(0)) ' 

It is clear that we need to estimate the functions Bi(6) and B(S) in order to 
estimate 

n(«, r) - n(«, r/2) = n(z p „ 1/2) + o (||n(u, r 

where the above equality is a direct consequence of Proposition [ 

Before we can estimate the integrals in (I4.12j) and (|4. 13[) we need to introduce 
some notation for integration on the unit sphere. We parametrise the unit sphere 
in M 2 according to 

0Bi(O) = {&(0); 0G(O,27r)}, 

where £i(</>) = (cos(<^>), sin(</>)). Inductively we define, for k > 2, the polar coordi- 
nates 

"01, 1P2, -.,ipk-i) = (sin(<fo._i)£ fc _i(0,V>i, ...,ipk-2),cos(ip k -i))- 

The unit sphere in R k is then defined by 

0Bi(O) = {e fe _i(0,Vi, -,^fc-2); (0,2tt), ^ 6 (0,7r)}, 

modulo a set of measure zero. 

With this parametrisation an area element on the unit sphere becomes 

(4-15) cL4 = det[%i ^ ... 1 ##i-# fc -2, 
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where £k-i is considered to be a column vector 
k x (fc — 1)— matrix in (|4.15[) is 
(4.16) 

#i)cos(0)P* = t 



Somewhat more explicitly the 



- sin(V>)PjLi 2 cos( 
cos(0)P,t"i 2 

-sin(^i)Pf- 2 2 



dn(0i) cos(0)P/- 2 2 



cos 

COS 



(^fc-3)cOS(0)P| =1 2 ^ fc _ 3 

(V'fc-a) sin(9i)P] I - 1 2 3 . #fc _3 cos(^_ 2 ) sinO)P/- 



COs(l/'fc-3)P,tlj #fe -3 
COs(Vfc- 3 )P,t 2 2 ^ )fe _3 



-sin(V>fc-3)sin(^ fc _ 2 ; 




COs(V>k-2) COS(0)P J 'L 1 3 
3 

1 



where we have used the notation 

p/r^n^sin^). 

We will denote the matrix in (|4.16p by M. By the anti-commutativity 
rows in the determinant function we have the identity 

det(M) = sin fe - 2 (V'fe-3)sin fe - 2 (V'fe-2)x 







cos(^ fc _ 2 )P*li 
cos(V>fc_ 2 )P,t 2 3 



cos(V>fc-2) 003(^-3) 
— sin(i/)fc_ 2 ) 



ivity of the 



(4.17) 





N 




u u 

-sin(^_3)sin(^_ 2 ) cos(^ fc _ 2 ) cos(^_ 3 ) 
-sin(^-2) 



= sin fc - 3 (^_ 3 ) sin fe - 2 (iP k _ 2 )det(N) 

where N((j), ipi, .., ipk-i) is the (k— 2) x [k— 3) — matrix satisfying sm(i})k-3) s'm(ipk-2)nij 
rriij for 1 < i < k — 2 and 1 < j < k — 3. Notice that N is independent of ipk~3 an d 

i>k-2- 

In order to estimate Bi we will use the identity in (|4.17p to write, with k = n, 
Bi{5) = X{ Ps >o}XidA dBl{0) = 

JdB 1 (0) 



(4.18) 



p27T p7T p7T 



Jo Jo Jo 



p2ir p7r pit 



pit pTT 



X{ Ps >o}X 2 l \det(M)\dipk-2dipk-3--d<j) = 
X{ Ps >o}xf |sin" -1 (V>n-3) sm n (-ip n - 2 )\ dip n -2diJjn-3 \det(N)\di/j n -4 



(4.19) 

= -2' 



Jo Jo Jo Uo Jo 
We will need some further simplifications 

Bi(S) = - / X{ PS >o}^dA aBl{0) 
iasi(o) 

X{ PS >o}xf \det{M)\dip k -idi:k-2---d(j) 



p2lT pTT /*71 

Jo Jo Jo 



>{0,ir/2y 



rir/2 rir/2 

/ / X{p s >0} X iS n ~ 1 ' n ('<Pn-2,lpn-l)dll)„-ld'lp n -2 

10 Jo 



= -2" 



/(0,tt/2)"- 



/ X{ PS >0}X 2 iS n 1 ' n {lp n -2,1pn-l)dlp n -ldlp n -2 

Jam 



\det(N)\dif> 
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/ 

J (0,7T 



(0,vr/2)2\A( A «) 



XfesXO^i S($ n -%, 1pn-l)dlpn-ldlp n -2 
h.i&V) + ^2,i(<5,/i), 



|det(JV)|d0 



where 



^-^"(Vn^^n-l) = |sin"- 1 (V„-2)sin™(V„-i)| 
and = i 7 ' _1 ((0, /j,) 2 ) where F is the stereographic projection 

'cos(^ n _ 3 ) cos(^„-2) 



(4.20) 



F(l/) n - 3 ,lf> n -2) 



% sm(ip n -a) ' sin(V'„-2)sin('0„_3)_ 

If /i is small then A(fi) ~ (7r/2 — /i,7r/2) 2 , the exact form of is unimportant 
as long as A(n) contains a small neighbourhood of the point (tt/2, 7r/2). We choose 
the particular form of in order to simplify some calculations further on (see 
equation (|4.24p ). 

We will estimate I\^{8,fj) and l2,i(S,fJ.) separately for \S\ small. Fix a fi > 
such that \5\ << fi << 1. The value of p, is not very important and can be chosen 
universal, depending only on n in particular p < cl in (14.32j) . 

To estimate l2,i(8,fj,) we notice that 

Vps = 2(^1X1,^2^2, -, 5„- 2 a; 2 I -2 J (1 ~ -%n)- 
By our choice of polar coordinates we have that when tp n -i £ (0, tt/2 — p.) then 

X n = COs(lp n -l) > C(J>- 

This means that the gradient of p$ is bounded from below by a constant times p 
on its zero level set. It is therefore very easy to estimate l2,i(8,p) by means of the 
co-area formula. 

By the co-area formula it follows that for t £ (0,1) and with the notation qs = 
T n ~ 2 5x 2 

\j t h,i{8,t5)\ 



'{«- 2 -<)/(<ls-fi<-i)=t} 
In particular 



V 



dA 



{(■- 



-*=)/(?«-«a£_i)=t} 



(4.21) 



|Ja, 4 (£,M)--fc,i(M)l <<?— • 

A* 



We need to work a little harder in order to estimate I\ (5, p). We begin to prove 
a simple lemma that will allow us to do some integrations explicitly module 0( la- 
terals. 

Lemma 4.1. Let 4>° , tp®, ipQ, ^ e fixed. Furthermore we let p > be a small 
constant and 1 < i < n. We use polar coordinates Xi((f>, "4>i,tp2i —jV'n— 2)- 
W^e also assume that 

n-2 

(4.22) ^S J x j (<l> ,4ll>l...,il>l^,n/2,K/2) 2 >0. 
3=1 
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Then there exist a constant c > such that 



(1-C/i) / ^W> ,^,...,V„-2) 2 

Ja(p) 



X (X{p <s >O}(0°>?,-,V'n-2) -X{po>0}(</ , °,V'?,->n-2)) 5" 1 '™# n _ 2 #n-l 



< 



^0 



S ? (X{p,,(x)>O}(<A0i\ -,0n-3,0n-2) ~ X{ Po >0} (0° , 0? , -,0n-2)) rfi„_idx„ 



<(l + C/l) / Xi(</>°,1pi, n _3, ^«-2) 2 
•* Mm) 

X (X{p 5 >O}(0°,V'?,-,V'n-3,V'n-2) -X{p o >O}(0°,V'?,-,V'n-2)) S^ 1 ' 71 <ty n -3<ty n -2i 

where 

S itj = |sin I (-0„- 3 )sin- : '(-0„-2)| , 

Xi(0,t/>1,02, ...,Vn-2) : 



and £/ie se£ A is i/ie stereographic projection of the two dimensional spherical area 

{x((j) ,1pi, ...,-0„_ 3 ,-0n-2); (</> n -3,^n-2) G (fl"/ — M> 7r /2) 2 } 

under the projection x — > x. 

Remark: Assumption \4.22]) is non-essential and only made for deGniteness and 
the result still holds if 

n-2 

SjXj^ , 0?, 2 °, 0£-4* */2, ^/2) 2 < 0. 

3=1 

Proof: It is trivial that 1 — c/i < sin(^„_3) < 1 and that 1 — c/i < sin(-0 n _ 3 ) < 1. 
Therefore 

(4.23) 1 - Qj/i < S 1 ' 3 < 1. 

Use the change of variables 

/ , , x V C0S(^„_ 3 ) COS(0 n _ 2 ) \ f~ N 

(0„_3, 0„-2j -> ^-t^ r, -r-r; r^-n \ = (^-1,^-2) 

in 

Xi(</> o ,0?,...,0„-3,0«-2) 2 

x (X{p«>O}(0°>01, -,i>n-3,1pn-2) ~ X{p o >0} (0° > 01 , •», V'n-2)) 5"^ 1 '"dV'„_ 3 d'0„_2 

(4.24) 

0?, 0n- 3 , 0™-2) 2 (X{ PS (S)>0} - X{p (£)>o}) S'™ _4 '™ _2 di n _i dx„ , 

it is in this change of variables that we use the rather awkward definition of A(fi) 
in order to get a nice area of integration to the right. 



Since yY^j=i x ] = sm (0n-3) si n (0n-2) we may estimate 
(4.25) (1 — cu.)xi < Xi < Xi 
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Notice that since 

n-2 

]T Sjx^ , v?, v 2 °, €-4, t/2) 2 > o. 

the integrand is non-negative so we may use (|4.23l) and (I4.25P in (I4.24[) to deduce 
the desired estimates. □ 

Lemma 4.2. Let \S\ << [i << 1. Also denote 

n-2 



i=i 

and 4>° , ip®, V-'„-4 fixed constants. Then, for i = 1, n — 2, 

ai? (X{ P5 >o}( ( / , >i ; --- I '0n-2) -X{ Po >o}) S n - 1 ' n d^ n - 3 dip n - 2 = 

I + 0(n) q s {<j) Q ,ifj^...,ir/2,Tr/2) , , 

~ |ln(|g 5 (0^V'i,---,7r/ 2 > 7r / 2 )l)|+O(|(5|//i+^|'5||ln(^|)|) 

* l — o 

and /or i = n — l,n we have 

A (X{p s >O}(0,V'lj-)^n-2) -X{p o >0}) <S n_1 ' fl #n-3*/>n-2 

(7r/2-p.7r/2) 2 

= 0(1*1) 

Proof: By Lemma 14.11 it is enough to prove the estimate for 
(4-26) / / x\ (x{ PS (x)} ~ X P0 (x)) dx n ^idx n , 



where J2"=i x ) = 1 - 

To simplify notation we will write 

K = qs(x). 

And we will assume that k > 0, if k = then the argument is simple and the case 
K < is treated analogously. 
Notice that 



X{P5(a:)>0} ~ \ else. 
For t=l,...,n-2we may write (|4.26|) as 

= 1-^ in( K )i 4 2 + 2(1 ^ 2) + o(|^|/m + mI*II Ml*l)l), 

where we have used the identity 

y \/T+ x 2 dx = ^z\A + + i In ^a; + V 1 + a; 2 ) 
to evaluate the integral. 
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For i = n — 1 we can calculate 







k + (1 + S)xl_ 1 - x n _ x dx n - Y = 0{n k) 



Finally, for i — n we get 



•^n (x{ps(x)} Xpo(x)) dx n —idx Tl 



yfn+(l-S)xl_ 1 



X n dx n 



x n dx n 



dx n -i = 0(k/j 2 ). 



a 



Proposition 4.3. If \5\ is small enough and C;(<5) is defined according to 

d(S) = Bi(S) - Bi(0) 
then there exists a universal constant c such that 



1 n-2 

~|<51n(|<5|)| < \ c i( 6 )\ < c|«5In(|*|)|. 



Moreover, if Si > Sj then Ci (S) < Cj (5) . 
Proof: In (|4.19l) we showed that we can write 

Bi(S) ~ Bi(0) = [hMl*) ~ ^,<(0,A*)] + MS,/*) - A(0,Ai)] 
We also showed, (|4.2ip . that 

[h t i{^^)-hM^)] = 0{\8\/n). 
Also in (14. 19)) we showed that we can write 
(4.27) Ji(*, M )-/i(0,m) 



(X{ps>0} - Xpo>o) S n 1 '"# n _ 3 #n- 



det(N)dA dBn -2 (0, ^„_4)- 



Furthermore we showed, in Lemmas 14.11 and 14.21 that the inner integral in (|4.27[) 
satisfies 



{X{ PS >o} ~ X Pn >o) S n 1 - n d%\) n - Z d^ n _- 1 

(l+0(p)) / xj (x{p s (x)>Q} ~ X{po>0}) dXn-xdXr, 



1 + 0(ji) 



qs(xi,..,x n -2) || ln(\q s (xi, ^ n _ 2 |)|| + + mI^I I m (l^l)l) 



for (x\, a; n _2) £ <9-B" 2 . Disregarding lower order terms we may conclude that 
(4-28) I M (<y,/i)-Ii,((0,/i) 
1 



g 4 1 ln( 9tf )|det(JV)(L4 aB „-2 + 0(|<J|/m + /#|| ln(|<5|)|). 



SB'; 



Let us denote the integrand F(qs), that is F(t) = t\ ln(|t|)|. We may estimate 
(4.29) \F(q s ) - \S\\]n(\S\)q s \\ < \SqsH\q s \)\ 
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where qs = mQs- Since qs is a second order polynomial with coefficients bounded 
by one it directly follows that 



(4.30) 



8B" 



\6q s hx(\q s \)\det(N)dA aB n-2 



= 0(1*1). 



By (jOU|) . (j4~29t and (|4T28|) we may estimate 
(4.31) hMri-hjfaiJ,) 

\H\s\) 



f)B" 



x dB 

We define the linear functional L : R ra_2 

J dB , : -2 q s x\dA dB n 

LS= ' 



q s det(N)xidA aB? -2 + 0(\5\/n + p\6\\]n(\6\)\). 



by 



Lb?- 2 qsxl_ 2 dA 



Writing L in matrix form we get 

L = XJ + X2J 

where Ai, A2 > 0, I is the identity matrix and 

1 1 1 ... 1 
1 1 1 ... 1 



J 



1 1 1 
1 1 1 



It is easy to see that u % = [1, 1, 1, 1] is an eigenvector corresponding to the 
eigenvalue Ai + (n — 2)\2 and that u J = e\ — ej for j = 2, n — 2 are eigenvectors 
corresponding to the eigenvalue Ai. In particular L have (n — 2)— linearly indepen- 
dent eigenvectors that correspond to strictly positive eigenvalues. We may conclude 
that det(L) > 0. It follows that there exist a universal constant cl > such that 

(4.32) \LS\>c L \5\. 

To finish the proof we notice that 

n-2 n-2 

£ iawi = E - ^(o)i = 1 h\s\)\\ls\ + o(\5\/» + m\H\s\)\) 



3=1 



3=1 



>-\6\\M\6\)\+0(\6\/v + »\8\\Hm)- 



And 



n-2 



n-2 



^2\C i {S)\=^2\B i (S)-B i (0)\ = \ H\S\)\\L6\ + 0(\5\/ll + n\5\\H\$\)\ 



3=1 



<c|<J||ln(|<J|)|+0(|(J|/^ + Ml*IIMI*l)l)- 
The proposition follows for /i small enough if \S\ << /i. 
The final statement follows easily since Ai > 0. 



□ 
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Remark: We will also use the notation C{§) = B{5) — 5(0). Notice that 

n 

(4.33) C{5) = Y,Ci{5) 

i-i 

since xf = 1 on the unit sphere. 



5. Proof of The Main Theorem. 



In this section we prove Theorem 11.21 
By assumption we have 

u(rjx) x n-l ~ X n 



lim 



rT-^o ||u(rja;)||i2( Bl ) Wx 2 ^ - x 2 || L 2 (Bl) 
for some sequence rj — > 0. Therefore 



(5.34) lim " ( "' rj "'°l, = 4-1 " 4- 



For any r > we can define a 5(r) according to 

II(u,r,Q) 



:pi( r )(£C). 



sup Bl |n(u,r,0)|| L 2 (Bl) 
With this notation (|5.34[) implies that (see I4.1l|) 

l%)l^0 

so we may, by choosing j large enough, assume that S(rj) is as small as we need. 
Also, from (|3.6j) and (|3.5j) we may deduce that 

sup |II(u, rj, 0)| > c|ln(rj)| 

Bi(0) 

for j large enough. 

If we denote sup B / ) \H(u, s, 0) = t s ss c| ln(s)| for s small enough and r 2 -j s is 
increasing in j (Lemma I3.3j) . Then Proposition 13.41 implies that 

(5.35) n(«, rj/2, 0) - II(u, r, , 0) + II(Z Pl , 1/2, 0) + 0(t- q ). 
The main step in our uniqueness proof for blow-up limits is 

Lemma 5.1. Let u be a solution to U.l\) and assume that ^"i^ — tttt = Pair\ 

1 v sup Bl(0) |n(u,r,0)| 1 0l y r ) 

for some 5(r) satisfying \S(r)\ < kq for some universal Kq. 
We also assume that 

n 

(5.36) Y,Ci(6[r)) < 0. 

i=l 

Then for each 7 < 1/8 there exist a constant C 7 such that if 

(5.37) max(*i(r),&(r) > ...,a B -2(r)) > C 7 r"T 
then 

max(^(r/2),fe(r/2),..,^ 2 ( r /2)) max (ft(r), <5»_ 2 (r)) 
l-5(r/2) 1-S(r) 
Moreover, if 5j < and 

8j < min (fc(r/2), &(r/2), *„_ 2 (r/2)) 
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then it follows that 



Sjjr/2) ^ SJr) 



(5.39) J \' ' < 

l-5(r/2) l-5(r) 

provided that \5.37\ ) holds. 

Remark: IfY^i=i Ci{o~(r)) > a similar result holds and the proof goes through 
with trivial changes. 

Proof: From (|5.35|) and (|4.14j) we can conclude that the coefhcient of the a^— term 
in n(u, r/2, 0) is 

(5.40) TrSjir) + K (n 2 Bj(5(r)) - nB(S)) + 0(r,r^). 
Next we make the following claim 

Claim: For j = 1, n - 2 we have n 2 Bj(0) - nB{0) = 0. 

Proof of the claim: This is easy to verify since we can calculate Z po , and thus 
Bi(0) explicitly (cf. [2j Lemma 4.4]): 
Define v : (0, +ooj x [0, +oo) ^ M by 



v(x n — % , x n ) 

7T 

2 V x n —i 



-4a;„_ia; n log(a;„_ 1 +x n ) + 2(x„_ 1 -a; n J --2arctan -nVi + I J' 



Moreover, let 

!w(i„_i,a; n ), x„_ix„ > 0,x„_i ^ 0, 
-^(-^n-i,^), X„_i < 0,X„ > 0, 
-u(x„_i,a; n ), s„_i > 0,a;„ < 0, 

and define 

- w(x n _i,a: ra ) - 7r(x|_i + ic|) + 8x ra _ia; n 

57T 

In particular, Z XnlXn (x n -\,x n ) is a rotation of Z po . It is clear that 



or equivalently 



U(Z Xn _ lXn , 1/2,0) = ^-tx n -ix ni 

IT 



U(Z Po ,l/2,0)= ] ^-(x 2 n _ 1 -x 2 n ). 



It follows that n Bj(0) — nB(0) = for j = 1, n — 2. This proves the claim. 

By the definition of Cj(8) we may thus write, for j = 1, n — 2, the coefficient 
of the Xj— term in II(u, r/2, 0) (that is equation (|5.40[) ) 

T r S,(r) - K (n 2 Cj(S(r)) - nC(S)) + 0{t~ 2 "I). 
Similarly we can express the x 2 l _ 1 coefficient of II(u,r/2,0) according to 

tv(1 - 5{r)) + ^ - K Q (n 2 C n ^(S(r)) - nC{8)) + 0{r^). 
The quotient of the x 2 and the x 2 l _ 1 coefficients of n(w, r/2, 0) is thus equal to 
TrSjjr) - Kp (n 2 C 3 (S(r)) - nC(S)) + Q{j-^) 
r P (l - 5(r) + ^-K (n 2 C n -i(6(r)) - nC(6)) + 0^) ' 
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Let us first prove the Lemma under the assumption 
(5.41) 5j(r) =max(*i(r),fe(r),...,*„_2(r)). 

Then the claim of the Lemma is 

T r 8 3 {r) - K (n 2 C 3 (S(r)) - nC{8)) + Q(r~^) ^ 8 3 (r) 

T r (l-~5(r)) + l -^-K Q (n2C n - 1 (6(r))-nC(S)) + 0(T- 2 '>) l-~6(r)' 
The inequality ([5^21 hold if 
(5.43) -K (l-8(r))n 2 Cj(8(r))+K n{l-S(r)-8j(r))C(S)+O (\8\Sj + r" 2 ^ > 0. 
From (15.41)) and Proposition 14.31 we have 

n— 2 n 

(n - 1)0,- (*) < Ci{8) + 0(\8\) = Ci(5) = C(S) 

i=l i=l 

where we used Lemma 1431 in the first equality and (|4.33[) in the last equality. Using 
this and Sj > in (|5.43l) we can deduce that the Lemma holds if 

-K (l - S)Cj(S) > O (\S\Sj + t- 2 i) , 

or equivalently if 

-Cj(5)>0{t-^), 
where we used that |Cj(5)| ~ \S\\ ln(|5|)|. 

In particular if \8\ is small and f|5.41[) holds then (|5 -38[) holds if Sj > C 1 t~ 1 . 
This is exactly what we wanted to prove. 

Next we chose any Sj < in order to prove (|5.39[) . 
Then the claim of the Lemma is 

r r S 3 (r) - Kp (n^jSjr)) - nC(5)) + 0(t£^) ^ Sjjr) 

T r (\-8{r)) + ^-K a (niC n - 1 (5{r))-nC{8)) + 0{T- 2 ' 1 ) 1-S(r)' 

The inequality (|5.44|) hold if 

(5.45) -K (l-8(r))n 2 Cj(S(r))+K n{l-S(r)-Sj(r))C(S)+O {\S\Sj + r"^) < 0. 
We either have that 

(5.46) C{8) < -CV" 7 
or 

(5.47) C 3 (S)<Cr^ 

for some universal C. This since if 8k — max (<5i(r/2), 82(r/2), 8,^2(^/2)) > 
C^t" 7 then Ck{8) < — cC 1 t~ 1 \ ln(r r ) so if C{8) > —Ct^ 1 then at least one of 
Cz(<5), for I = 1, n — 2, must satisfy Ci{8) > cC 1 t~ 1 \ In(r r )| >> C^t" 7 since 
\8\ « 1. By the monotonicity of Ci(8) it follows that Cj(8) > Ct' 1 . 

In either case (|5.46l) or (|5.47[) it follows that (|5.45[) holds true. The Lemma 
follows. □ 

We may now proceed with our proof of the main Theorem. From Lemma 15.11 
and (11.21) it follows that 



(5.48) \8{r)\ <Ct-i. 

If not then we have by Lemma 15.11 that 

rnax(<5 1 (r/2),5 2 (r/2),..,5 n _ 2 (r/2))>m aX (5 1 (r) ) 5 2 (r),..,5 n _ 2 (r)) 



17 



if 

max{Si(r),S 2 (r),...,8 n ^ 2 (r)) >0 

and 

miii(<Si(r/2),(S 2 (r/2),...,(S n _ 2 (r/2)) < min (<Si(r), 5 2 (r), ...,(S„_ 2 (r)) 

if 

min((5i(r),(5 2 (r),...,(5„_ 2 (r)) < 0. 

Since T r / 2 k > T r i 2 i for k > I we may iterate this and conclude that if (|5.48l) is not 
true then 

lim max {5 x {r/2% 5 2 (r/2 k ), <5„_ 2 (r/2 fe )) > max (5 x (r), S 2 (r), <5 n _ 2 (r)) 
and / or 

Urn whx(5i{r/2 k ),5 2 {r/2 k ),...,5 n - 2 {r/2 k )) < min(&(r),&(r), J„_ 2 (r)) . 

fe— >-oo 

This would contradict ()1 .2|) . 

So (|5.48[) has to hold. This implies in particular that 

n(u,r,0) n(«,r/2,0) 
sup Bl |n(u, r, 0) | sup Bl \U(u, r/2, 0) 

We may iterate and conclude that 

II(u,r,0) n(u,r/2 fc ,0) 
sup Bi |n(u,r,0)| " sup Bi \U(u,r /2 k ,0) 



< C r- < Ct- 1 - 1 . 



<cJ2cr;; 2 ^< 



(5.49) <C^(feln(2)+ln(l/r)) 



-1-7 



since r r > c|ln(r)|. Since 7 > it follows that (|5.49p is convergent and we may 
directly conclude that 

u(rjx) 
hm — 

r-»o \\u(rjX)\\ L 2 {B ) 



exists. The first claim (|1.3|) of Theorem II .21 follows. 
That 



SnB ro (o)nL; |>? <7?«+4)| 



consists of two C 1 manifolds intersection at right angles at the origin is now standard 
(see Corollary 9.2 or in [3]). 
To prove that 

5„_ 2 ns ro (o) 

is contained in a C 1 manifold of dimension (n — 2) for some small ro we may proceed 
as in Theorem 12.2 in (3J. This proves Theorem 11.21 
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